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The merger of two spiral galaxies is believed to be one of 
the main channels for the production of elliptical and early-type 
galaxies. In the process, the system becomes an (ultra) lumi¬ 
nous infrared galaxy, or (U)LIRG, that morphs to a quasar, to a 
K+A galaxy, and finally to an early-type galaxy. The time scales 
for this metamorphosis are only loosely constrained by observa¬ 
tions. In particular, the K+A phase should follow immediately 
after the QSO phase during which the dust and gas remaining 
from the (U)LIRG phase are expelled by the AGN. An inter¬ 
mediate class of QSOs with K+A spectral signatures, the post- 
starburst QSOs or PSQ, may represent the transitional phase be¬ 
tween QSOs and K+As. We have compiled a sample of 72 bona 
fide z < 0.5 PSQ from the SDSS DR7 QSO catalogue. We find 
the intermediate age populations in this sample to be on average 
significantly weaker and metal poorer than their putative descen¬ 
dants, the K+A galaxies. The typical spectral energy distribution 
of PSQ is well fitted by three components: starlight; an obscured 
power-law; and a hot dust component required to reproduce the 
mid-IR fluxes. From the slope and bolometric luminosity of the 
power-law component we estimate typical masses and accretion 
rates of the AGN, but we find little evidence of powerful radio- 
loud or strong X-ray emitters in our sample. This may indicate 
that the power-law component originates in a nuclear starburst 
rather than in an AGN, as expected if the bulk of their young 
stars are still being formed, or that the AGN is still heavily en¬ 
shrouded in dust and gas. We find that both alternatives are prob¬ 
lematic and that more and better optical. X-ray, and mm-wave 
observations are needed to elucidate the evolutionary history of 
PSQ. 

Galaxies: starburst; QSOs. 


1. Introduction 

The possibility that at least some elliptical galaxies could be the 
end result of the merger of two spirals of similar masses was 

Send offprint requests to: Jorge Melnick, e-mail: jmelnick@eso. org 


first proposed by Toomre & Toomre ( \912) (see also Toomre 


1977 1 and has since been explored in numerous investigations 
(see e.g. [Hopkins et al. 2008a for a recent review). While in¬ 
triguing, the idea of the metamorphosis of spirals into ellipti¬ 
cals has been questioned on dynamical and photometric grounds. 
Nevertheless, one by one these objections have been resolved 
by increasingly sophisticated observations and numerical mod¬ 
elling. 

The main difficulty lies in suppressing star formation in the 
merger of two gas-rich systems to return a quiescent remnant. 
This process (quenching) may be aided by feedback processes, 
and especially by the activation of an active nucleus (expected 
if mergers drive gas and dust to the centre by tidal torques 
(e.g. Mihos & Hcrnquist; 1996[ >. In the detailed simulations by 


Hopkins et al. (2008a|bl the merger of two gas-rich systems pro¬ 


duces an AGN that eventually stops or prevents further star for¬ 
mation. The resulting remnant goes through a brief AGN phase 
and is then identifiable as a post-starburst (K+A) galaxy as the 
initial star formation episode fades away. A recent investigation 
has revealed that a median 50% of the stellar mass in a sample of 
808 K+A galaxies has intermediate ages and high metallicities 
(Melnick & De Propris 2014| >, consistent in fact with the metal¬ 
licities of elliptical galaxies of similar masses. Thus, the present 
day stellar populations of at least some elliptical galaxies appear 
to have been formed in a colosal burst of star formation induced 
by the merger of two gas-rich galaxies of similar masses (the 
preferred model for the origin of K+A galaxies Snyder et al. 

|20TTT >. 

The signatures of this process may be difficult to identify. 
There are several lines of evidence now favouring a secular sce¬ 
nario for AGN evolution (e.g. Jahnke & Macci6|[2011] l. X-ray 
selected AGN hosts in deep HST imaging are not observed to 
be more disturbed than a similar sample of inactive galaxies at 
the same redshifts (e .g. |Cistemas et aLl|2011[ |Kocevski et al.| 
201 2[ Villforth et al.||2014|and references therein), while both 


De Propris et al. ( 2014[> and Scott & Kaviraj (2014 1 do not find 


evidence for a strong enhancement of the merging fraction in 
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closely interacting pairs, despite the presence of tidal signatures. 
Similarly, Teng et al. ( |2012[ i searched for AGN activity in com¬ 
panion galaxies to QSO hosts and detected just one AGN in their 
sample of 12 galaxies. 

It is of course possible that early QSOs are highly obscured: 
Koss et al. ( 2011 1 > shows that hard X-ray selected AGNs, which 
are often Compton thick, are hosted by highly disturbed galax¬ 
ies, while Urrutia et al. (2008) find that a large fraction of Type II 
(red, dust obscured) QSOs lie in morphologically disturbed sys¬ 
tems. Canalizo & Stockton ( 2013} ) argue that QSOs with strong 
absorption lines have evidence of tidal interaction, but this is 
only present at low surface brightness levels, which argues for 
gas-poor, low mass, or even old remnants, that may or may not 
be connected with the observed present-day activity. 

If we accept the QSO-merger model, we expect that at least 
some young QSOs may be associated with a post-starburst signa¬ 
ture as they have just suppressed star formation and are emerg¬ 
ing from their dust-enshrouded phase. In a previous paper (De 
Propris & Melnick] |2014| ) we have shown that no K+A galaxy 
at any age hosts an AGN, although there is clear evidence of 
feedback (see also Wong et al. [2012 for a similar conclusion). 
Therefore, suppression of star formation and QSO activity must 
rapidly follow each other (see Teng et al .|2012 1. Thus, the par¬ 
ticular class of quasars known as post-starburst quasars (PSQ; 
see [Cales et al.|20r3| for a recent review), would be the immedi¬ 
ate ancestors of K+A galaxies and therefore the descendants of 
recent mergers. 

Observationally, PSQ are characterised by strong emission 


not sufficient to ascertain the presence of strong H6 in absorp¬ 
tion. Finally we ran the Starlight code (Cid Fernandes et al. 
2005) to model the stellar populations of these 83 objects using 
exactly the same prescriptions of Mel nick & De Propris | ( |2014| ) 
for the K+A galaxies: BC03 (Bruzual & Charlotf 2003} stellar 
library for 25 ages and 5 metallicities. From the Starlight fits 
we found that the SDSS spectra of 10 out of our 83 bona fide 
PSQs were of insufficient quality to derive reliable stellar popu¬ 
lations, either because the S/N was too low (S /N <2 at 4100A), 
or due to significant gaps at some critical wavelengths. 

The essential statistical properties of our final sample of 72 
PSQ are presented in Figure [Tj in comparison the properties of 
the K+A galaxies from Melnick & De Propris| ( j2013] ). The com¬ 
parison is relevant because the main purpose of our investigation 
is to establish whether PSQ can be the progenitors of K+As. 

Figure [T] shows that - as QSOs - PSQ are rare in the local 
Universe, even compared to K+A galaxies. On the other hand, 
the stellar masses of PSQ appear to be comparable to, but some¬ 
what larger than, K+As. We note, however, that the PSQ masses 
are rather uncertain given the generally low S/N of our spectra. 
More surprising is that the Petrossian radii of the two samples 
are comparable given the differences in redshift. While this is 
probably due to the relatively low resolution of the SDSS im¬ 
ages, it confirms that, as discussed by Melnick & De Propris 
(2013]), the photometry and spectroscopy of our PSQ sample 


lines and strong absorption of the high Balmer lines. Cales et al. 
( 2013| > concluded that PSQ with early-type hosts are likely the 
result of major mergers, so this class of objects may offer a way 
to witness the final stages of the metamorphosis of spirals into 
early-type galaxies. In this paper we present a systemic study of 
the properties of a sample of 72 PSQ from the SDSS selected 
from the DR7 Catalog of Quasars (Schnei der et al.| 2010} ), par¬ 
alleling our study of K+A galaxies. We use the SDSS spectra 
to derive the stellar populations, which are then used to predict 
the broad-band SEDs from the far ultraviolet (FUV) to the mid 
infrared (mid-IR(. Throughout this work we adopt the following 
cosmological parameters: Cl M = 0.3, O a = 0.7, and // () = 70km 
s 1 Mpc -1 . 


2. Data and models 

The DR7 Quasar Catalog (DR7Q) contains data for 105,783 
SDSS quasars selected as described by Schneider et al. ( 2010} >. 
We downloaded the SDSS spectra of the 8492 objects in the 
DR7Q with z < 0.5, which we searched for signatures of post- 
starburst stellar populations using a slightly more relaxed crite¬ 
rion than that used by Goto ( 2006) 1: we looked for objects with 
(rest frame) equivalent withs of H<5 greater than EW(Hd) > 3A 
in absorption. We found that about 1.5% of the objects (125 out 
of 8492) satisfied this criterion. 

Our PSQ fraction is significantly lower than the fraction of 
4.2% PSQs found by Goto] ( [2006} in a volume limited sample of 
SDSS galaxies. This is not surprising given that Goto gives the 
fraction of all AGN that have strong Il y in absorption, while our 
sample is restricted to objects classified as QSOs, that is, only 
to AGN with broad emission lines. In fact, the three examples 
shown by Goto ( 2006| > are classified as (Seyfert 2) galaxies in 
the SDSS and therefore are not included in the DR7Q. 

We visually inspected the spectra for spurious detections, 
which led to a further reduction of the sample to 83 objects; 
while some were clearly not PSQs, the S/N ratios of others was 


are not significantly affected by systematic aperture matching 
effects. 

The star formation histories for our 72 bona fide PSQs were 
combined to build the spectral energy distributions (SEDs) of 
the stellar components from the FUV to the mid-IR using the 
BC03 stellar library to convert masses and ages to photometric 
fluxes, which we compared to the photometric observations. We 
emphasise that we do not use the spectral fluxes to construct the 
SEDs, but the stellar populations derived from the spectra. This 
allows us to essentially extend the spectral coverage to the FUV 
and mid-IR. 

The essential characteristics of the broad-band photomet¬ 
ric dataset are: optical (ugriz) from the SDSS (Abazajian et 


al. 2009 1 ; Vacuum UV data are from the GALEX imaging 


databases (Morrissey et al.| |2007|); Near-infrared data in JHK 


come from 2MASS (Skrutskie et al. 


data from the WISE (Wright et al. 


2006)); and Mid-infrared 


2009 1 datasets. The pho¬ 


tometry and the input SDSS spectra were corrected for Galactic 
foreground extinction using the values of given in the DR7Q 


(Schneider et al. 2010) using the standard Galactic reddening 


law ( |Cardelli et al.|1989) henceforth CCM). 

The synthetic spectral energy distributions are calculated as 


25 5 


m )= 


(i) 


«=I J= 1 


where M lf are the BC03 models for the i = 1,25 age bins 
and j = 1,..,5 metallicity bins, and S tJ are the mass fractions 
returned by Starlight for the relevant age and metallicity bins. 
The photometric zero-point k is independent of wavelength. 


3. Results 

3.1. Stellar Populations 

Figure |2]presents in the left panel the average stellar populations 
of the full sample of 72 PSQ analysed with Starlight. These 
plots give fractions of the total stellar mass as a function of age 
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Redshift 



Fig. 1 . Statistical properties of the PSQ sample (in black) compared to those of their putative descendants the K+A galaxies (in red). 
The number of objects in each bin has been divided by ten for the K+A radii and by five for the masses.’ 


for three broad metallicity bins (metal poor; solar; over-solar) 
chosen to make the plot less cluttered. For comparison, the right 
panel presents the mean stellar populations of our sample of 808 
K+A galaxies from Melnick & De Propris (2014). Two remark¬ 
able differences are immediately obvious: (1) the average frac¬ 
tion of mass in intermediate age stars is about 15% in PSQ com¬ 
pared to ~ 45% for K+A galaxies; and (2), while in K+As the 
vast majority of the intermediate-age stars are metal rich, in PSQ 
this is the case for only a small fraction of the stars. 

It is notoriously difficult to estimate formal statistical er¬ 
rors for the stellar populations derived from population synthesis 
modelling (see e.g. Melnick & De Propris|2014 for a discussion 
in the case of K+A galaxies and Starlight). Systematics domi¬ 
nate over formal fitting errors, in particular the well known de¬ 
generacies between age and metallicity. In the case of PSQ one 
should be particularly concerned by the effect of a featureless 
power-law continuum from the AGN. 

Starlight offers various possibilities to include a power-law 
continuum, and in particular, to deal with different extinctions 
for different populations. We ran many experiments to establish 
whether a combination of power-law and extinction could affect 
the resulting stellar populations. We found that, while including 
these combinations may change the detailed distribution of ages 
and metallicities, the total fractions of intermediate age stars re¬ 
mained substantially unchanged. 

We obtained very consistent results by either including a 
power-law component of slope a = -1.5 (/) (/t/0.55/r)~ L5 ) 

reddened by Ay — 1.3, as determined from our SED fitting 
discussed in the next section, or by including an un-reddened 
power-law of the same slope and allowing Starlight to choose 
the best fitting extinction for each object. The former is shown 
in Figure [2] and is the one we used to compute the total stellar 
masses, but the difference between the two methods are surpris¬ 
ingly small. 

From these and the other experiments we ran, we are confi¬ 
dent that the differences between the stellar populations of PSQ 
and K+As shown in Figure[2]are not due to our choice of how to 
model the power-law continuum. 

Another systematic effect is contamination of the Balmer 
stellar lines by nebular emission, which could result in lower 
fractions of intermediate-age stars. Visual inspection of the spec¬ 
tra shows this effect to be present in some of our objects, but even 
at Hb (we mask H a,j3, y) the contamination is very small. In or¬ 
der to quantify this effect we subtracted the population synthesis 


model from the original spectra to generate a pure nebular spec¬ 
trum. We then subtracted the smoothed nebular spectrum from 
the observed spectrum to generate an emission-line free spec¬ 
trum, which we fitted again with Starlight. 

This procedure inevitably adds noise to our already rather 
noisy SDSS spectra. We have tested the procedure with rather 
discouraging results: we get no net increase in the fraction of 
intermediate-mass stars and the metallicities do not change sig¬ 
nificantly. Masking out Hh in the original spectra gave similar 
negative results. We also tested stacking the spectra of all the 
objects that allows us to remove the emission lines from the 
stacked spectrum without adding noise. The change in the stellar 
populations with and without removing the emission lines was 
negligible and EW(Hb) changed by less than 10%. 

We performed a number of other tests to evaluate the sys¬ 
tematical errors in our results: changing spectral range to elim¬ 
inate the noisiest spectral regions; changing the reddening law; 
and masking various spectral features. While the stellar popu¬ 
lations change in detail, the overall fractions and metallicities 
of the intermediate-age component never changed significantly 
and in particular never increased. We are confident, therefore, 
that the discrepancy between PSQ and K+As is not due to errors 
in the fitting procedure. 

3.2. Synthetic SEDs 

Figure [3] compares the average synthetic SED computed using 
Equation 1 for the mean of the Starlight stellar populations 
(shown in Figure [2| to the observations. In order to remove the 
intrinsic variance in luminosity we normalised the fluxes to the 
SDSS /'-band ( Fansium )■ The observed fluxes show a clear ex¬ 
cess above A ~ 0.5/r, which is very different from the case of 
K+As where the excess is only observed in the mid-IR bands at 
A > 5/j. In order to fit the PSQ, therefore, in addition to a hot-dust 
component parametrised as a black-body of temperature T c i ust as 
we did for K+A galaxies, we added a reddened power-law com¬ 
ponent as expected for AGN. 

We emphasise that the error bars plotted in the figure are not 
observational photometric errors, but the rms scatter of the nor¬ 
malised fluxes about the mean of the sample. The mean obser¬ 
vational errors range from 0.26 magnitudes at 22jj and 0.23 mag 
in FUV to < 0.03 mag in the SDSS g' to z! bands. Thus, while 
PSQ form a remarkably homogeneous class of objects at optical 
wavelengths, this is not the case in the UV and IR bands. This is 
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Fig. 2. Comparison between the stellar populations of PSQ on the left and of K+A galaxies on the right for three broad metallicity 
bins. 


different from K+A galaxies that have remarkably similar SEDs 
in the optical and in the IR up to 4.5p (jMelnick & De Propris 
12014 ). 

We thus need six free parameters to fit the observed SEDs: 
Tciust', (y uv - the slope of the power-law component; and ,4[’/' - 
the visual extinction of the power-law, plus the photometric zero- 
points of the three components: stellar populations; obscured 
power-law; and hot-dust. An additional free parameter is the red¬ 
dening law for the stellar and power-law components, which is 
very relevant for the UV bands. 

In principle the extinction is a free parameter only for the 
power-law component. In the case of K+A galaxies we corrected 
the photometric observations using the extinction values from 
the Starlight fits. However, this is not possible, or at least is not 
straightforward, in the case of PSQ because the power-law com¬ 
ponent contributes a significant fraction of the continuum flux 
in the SDSS spectral band. We verified that this does not change 
the stellar populations significantly, but it certainly affects the re¬ 
sulting extinction, which in turn affect the SED fits and thus the 
optical power-law continuum. We tested an iterative procedure, 
which is very time consuming and does not easily converge. 

We therefore followed a different procedure to estimate the 
extinction corrections for the stellar component, 4 *. First we 
fixed the photometric zero points of the power-law and hot-dust 
components by forcing each to fit the WISE 4.5p and 12/; bands 
respectively; the power-law completely dominates the 4.5p flux 
and the Black-Body completely dominates the 12/j emission, 
with negligible extinction corrections. For the stellar component 
we fixed the zero point in the SDSS /-band where stars still con¬ 
tribute a significant fraction of the total flux while, not being 
negligible like in the mid-IR, the extinction corrections are still 
much lower than in the bluer bands. 

After fixing the photometric zero points we fitted remaining 
two free parameters (we fixed the dust temperature to T<i us , — 
250K ab initio ): the slope (ffyy) and the extinction A (7‘ of the 
obscured power-law component, and we computed the stellar 
extinction A* v by forcing the FUV flux to match the observed 
values (even for a small amount of extinction the contribution of 
the obscured power-law is negligible in the FUV). We then re¬ 
computed the zero points on the extinction-corrected fluxes and 


iterated until the zero point of the stellar component converged 
to better than 1%, which takes only a few iterations. 

Following our work in Paper I, we experimented with two 
extreme extinction laws: the canonical Galactic law of Cardeflil 
eFakl (fl989l> - CCM , and the UV-strong SMC/30Dor law of 
Gordon et al. ( {2003 i, which we refer to as GDI following the 


nomenclature of Starlight. Both laws are reasonably similar 
in the optical and near-UV, but differ significantly in the FUV. 
Also, while CCM shows the strong 2200A bump characteristic 
of the Galaxy, the bump is absent in GDI. We verified that GDI 
provides a much better fit to our PSQ data although, within the 
redshift range of our sample, the 2200A bump does not affect 
the observations. The photometric data uncorrected for extinc¬ 
tion are shown as blue squares in Figure [3] while the extinction- 
corrected data are shown in red, for the best fit mean stellar ex¬ 
tinction of Ay = 0.27 and the GDI reddening law. 

Figure [3] shows that our simple three component model pro¬ 
vides a remarkably good fit to the observations over more than 
two decades in wavelength from the FUV to the mid-IR. There is 
clearly some degeneracy in the parameters, notably between in 
auv and 4(7', although we find solutions only for -1.2 > auv > 
-1.8 and Ay L > 0.8. From a visual inspection of the+ 2 contours 
the best fit is obtained for auv = -1.5 ±0.2 and Ay L — 1.3 + 0.3. 
Interestingly, the power-law slope is constrained to values re¬ 
markably similar to the average value observed for X-ray lumi¬ 


nous AGN (f A oc /T , Wu et al. 2012). 


This reddened power-law component is the one we in¬ 
cluded in the Starlight population synthesis models described 
in Section l3~n 


4. Emission-lines 

4.1. Gas excitation 


Figure |4| shows the WHAN diagram described by Cid Fernandes] 


et al. ( [2010 2011[ ) as a powerful emission-line diagnostic 


method for objects with weak emission-lines, or systems, such as 
ours, where the higher Balmer lines are severely affected by un¬ 
derlying stellar absorption making diagnostics such as [OITI|/H/J 
(or [OII]/[OIII] that requires extinction corrections) rather un¬ 
certain. Since most of the objects in our sample have broad Her 
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Fig. 3. Best fit model SED for the average of the 72 K+A galaxies in our sample compared to the average observed broad-band 
fluxes corrected for extinction as described in the text and normalised to the SDSS /'-band The blue symbols show the 

observed photometric data points and the red symbols show the observations corrected for extinction using the visual extinction 
shown in parenthesis in the figure legend, and the GDI extinction law. The three components of the fit are shown colour-coded as 
follows: cyan shows the contribution of the stars computed from the Starlight population-synthesis fits to the SDSS spectra; green 
shows the obscured power-law; and the magenta represents the Tj ust = 250 K hot dust component, which is required by the mid-IR 
data. The sum of these three components is shown in black. The error bars are not observational errors but correspond to the lcr 
dispersions of the normalised fluxes for all the objects included in each band, the number of which varies from 45 objects for the 
FUV band to 72 for the SDSS bands. The reddened power-law component was included in the population synthesis models, thus 
the label Pop=R15 in the figure legend. 


emission lines, we used the package PAN in IDL (Dimeo, [2005| > 
to simultaneously fit the narrow (Hu and [Nil]) and the broad 
components. This was only possible for a subset of all our ob¬ 
jects because at their redshifts the SDSS spectra around Her are 
rather noisy due to contamination from atmospheric features, 
and in some cases outside the spectral range. 

As expected (the objects were selected from a QSO cata¬ 
logue) most of the objects in our sample for which we could 
measure the parameters fall in the region of strong (and presum¬ 
ably real) AGN (sAGN). A few are weak AGN (wAGN formerly 
LINERS), while only one PSQ appears to be powered by young 
stars (SF). 


4.2. Gas kinematics 

In the process of separating the broad and narrow components 
of Hu we noticed that for the majority of our objects the broad 
component appeared to be systematically shifted relative to the 
narrow components, mostly, but not in all cases, to the blue. This 
is illustrated in the left panel of Figure |5]that plots the radial ve¬ 
locity difference between the centroids of the broad and narrow 
lines measured from the multiple Gaussian fits. 

The right-panel of the figure shows the velocity shift be¬ 
tween the narrow components (in most cases measured from the 
two [OIII]4959,5007 lines) and the stellar lines as measured by 
Starlight from cross-correlation with the best fitting stellar tem¬ 
plates. Negative values correspond to gas moving towards the 
observer and positive values to gas moving away from the ob¬ 
server relative to the stars. From repeated fits and from the two 
[OIH] emission lines we estimate the (star-gas) radial velocity 


differences to be accurate to ±20km/s, so the figure shows a clear 
excess of blue shifted gas emission. 

In general the H/J emission line in PSQ is either completely 
absorbed by the underlying absorption, or broad and rather 
weak, but for some cases H/3 is strong and narrow showing clear 
P-Cygni profiles. Some of these cases are shown Figure [6] that 
show blue-shifted emission relative to the stellar line, confirm¬ 
ing the predominance of blue shifts for the emission lines. 


4.3. Black-Hole masses and X-ray luminosity 


The simplest way to determine Black-Hole masses for our PSQs 
is either from the Eddington luminosities, or from the masses 
of the spheroidal component of the host galaxies. The statisti¬ 
cal study of X-ray luminous AGN by Wu et al.|( 2012| provides 
a useful way to estimate the Eddington ratio r\ E dd = /-«„// L E dd 
from the slope of the power-law continuum in the UV. The av¬ 
erage power-law slope for our sample auv = -1.5 (in f t units) 
is typical for the QSOs in the sample of Wu et al. ( j2012| l and 
indicates an Eddington ratio r/ E dd ~ 0.3. 

Runnoe et al. (2012|> provide bolometric luminosity correc¬ 


tions for quasars at 0.145/r and 0.3/r. Assuming a power-law 
slope auv = -1.5 for all the objects in our sample, we can use 
the dereddened power-law luminosities at these wavelengths to 
estimate the bolometric luminosities using the recommended re¬ 
lations of Runnoe et al. ( 2012) >. We then use the Eddington lumi¬ 
nosity, LEdd = 3 x 10 4 Mbh in solar units, to compute Black-Hole 
masses M™ as 


M E B h = 


LboI 


3 X l0 4 T]Edd 


M Q , 


( 2 ) 
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Fig. 4. WHAN diagram for the 72 PSQ in our sample. The grey points show our sample of 808 K+A galaxies for comparison. 




Fig. 5. Left Radial velocity differences between the broad and narrow components of the Her line measured fitting multiple 
Gaussians. Right. Histogram of the difference in radial velocity between the emission lines (V em and the absorption lines (V a b s )- 
The measurement errors are about +20km/s. 


where as before rjEdd is the Eddington ratio. The resulting distri¬ 
bution of black-hole masses is shown by the black histogram in 
Figure [7] 

The second histogram in red shows the BH masses estimated 
from the mass of the spheroidal component M sph using the to¬ 
tal stellar masses M, from the population synthesis models as a 
proxy for M sp i,. For simplicity we assumed M sp h = 0.85M* for 
all our objects, and used the correlation between M sph and M BH 
from f 


Scott et al. 


to infer BH masses. 


(2013): logM* BH = 0.97log(M sph /3xl0 n )+9.27 


Two distributions in Figure [7] are in remarkably good agree¬ 
ment considering the approximations that we have made, in par¬ 
ticular using the same power-law and the same stellar mass to 
bulge mass ratios for all our objects. Still, the Eddington ratio 
r/Edd = 0.3 depends only weakly on the UV power-law slope (Wu| 


et al.||2012j), while if anything the assumption of M sp h = 0.85 A/., 
is rather conservative. 


The largest uncertainty arises because a substantial fraction 
of the total luminosity at 0.55 n is contributed by the power-law 
component (Fig. [3j. So, in order to calculate the total stellar 
masses of each object, we have assumed that for all the objects 
in our sample the power-law component contributes 55% of the 
visible light. This correspond to the ratio of stellar to power-law 
luminosity of L^/Lpi = 0.81 (Figure [3j. 


Cales et al. (2013) studied a sample of 38 PSQ with substan¬ 


tially better S/N spectra, and using the more standard method of 
combining line widths and optical luminosities. Their masses for 
12 objects in common are comparable, within the uncertainties, 
with our M* bh values, indicating that the assumptions that went 
into our calculations are not entirely unreasonable. 

An independent test is to use the X-ray luminosities de¬ 
rived from the bolometric corrections. Runnoe et al.[ (20T2| give 
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Fig. 6. Examples of some of the few objects in our sample that have narrow and relatively strong H fi emission lines. The P-Cygni 
profiles show that the gas is blue shifted relative to the stellar component. 


two separate relations for radio-loud and radio-quiet objects. All 
the PSQ in our sample are in the NVSS ( Condon et al.| 1998| > 
and, as expected, the vast majority are radio quiet. Only six ob¬ 
jects in our sample are radio loud with fluxes well in excess of 
3mJy/beam, which seems to be consistent with the fact that at 
similar luminosities only about 10% of all QSOs are radio loud 
(Padovani |1993 ). Inverting the bolometric correction for radio¬ 
quiet QSOs: logL Bot = 33.08 +0.291ogL 2 -iOA'cV we get the distri¬ 
bution of 2-10keV X-ray luminosities shown in Figure [8] 

At these luminosities we would expect only three of our 
objects to be detected by ROSAT, while none appears in 
the ROSAT all sky catalogue ( |Voges et al.[ \1999) , which at 
the mean redshift of our sample has a 3cr detection limit 
of logL 2 -\okeV ~ 44.3 erg/s. Nine of our PSQ are in the 
XMM or Chandra archives, but only two (J025938.15+004206.4 
and J074524.97+375436.6) are robust detections, and a third 
(J094820.38+582526.4) is an upper limit. The mean X-ray lu¬ 
minosity of these objects, logLx = 43.2 erg/s, is consistent 
on average with our estimates from the bolometric luminosities 
(logLx = 43.1). We can confidently conclude, therefore, that the 
AGN in PSQ appear to be normal from the point of view of their 
black-hole masses, and their radio and X-ray luminosities. 


5. Discussion 

5.1. Kinematics 

In the scenario of |Hopkins et al.| ( [2008a) i, the QSO phase is pre¬ 
ceded by the Antennae/ULIRG phases characterised by the for¬ 
mation of numerous starburst clusters and super-clusters in the 
arms of both spiral components. While large amounts of molec¬ 


ular material may be channelled to the nuclear regions by tidal 
forces during the merger, nuclear starbursts in these phases are 
either still weak, or heavily enshrouded in dust such that their 
emission is not visible at optical and near-IR wavelengths. By 
the time massive starburst activity is finally ignited in the nu¬ 
cleus, all star formation in the arms would have ceased leaving 
behind a large intermediate age population of stars distributed 
across the remains of the tidally distorted spiral arms. 

Therefore, while the emission lines in PSQ mostly arise from 
the nuclear regions (e.g. Sanmartim et a l.|2014) , the stellar lines 
originate mostly from the remnants of massive starburst clusters 
in the spiral arms, so it should come as no surprise to observe a 
wide distribution of radial velocity shifts between stars and gas. 
And, since these shifts depend mostly on the - in principle ran¬ 
dom - orientation of the merger with respect to the line of sight, 
we would expect to observe a uniform distribution of velocity- 
shifts centred around zero. However, as shown in the right panel 
of Figure [5] while we do observe a uniform distribution, the gas 
appears to be systematically blue shifted relative to the stars. 

On the other hand, we observe a rather messy distribution 
of relative velocities between the narrow and broad components 
of Ha, with the broad components either blue-shifted or red- 
shifted relative to the narrow lines. Therefore, if the broad-line 
emitting region (BFR) is at rest relative to the systemic velocity 
of the galaxy, the centroid shifts we observe must be due to a 
combination of extinction and inflow, outflow, or rotation of the 
BFR. 

For a spherical geometry, the predominance of blue shifts in 
our sample would indicate a net outflow of the BFR, contrary to 
the model of Gaskell & Goosmann ( 2013) 1, but other geometries, 
such as an inclined rotating disk plus an oblique obscuring torus. 
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Fig. 7. left. Black-Hole masses calculated with two different 
methods shown in the lagend. The histogram in black (lleft) was 
calculated using the Eddington ratios and the bolometric lumi¬ 
nosities. The histogram in red (right) shows the distribution of 
masses calculated using the relation between Mbh and M sp h. 
right. Correlation between and M* BH . The solid line shows 
a least-squares fit with coefficients indicated in the legend. 



log 1^(2—lOkeV) [erg/s] 


Fig. 8. Distribution of X-ray luminosities for PSQ calculated us¬ 
ing the bolometric corrections for radio-quiet QSOs shown in 
the legend. 


could also provide the variety of blue and red shifts observed in 
PSQ (Stern et al., 2015} . Beyond asserting that there are ways 
to explain the observations of PSQ, however, the kinematics of 
the BLR is a very complex issue that is well beyond the scope of 
this paper. 

Nevertheless, we may safely conclude that the distribution 
of the relative velocities between the narrow emission lines and 
the stars result from a combination of rotation plus a net outflow 


(in most cases) of the narrow emission line emitting gas. This 
conclusion agrees with the spatially resolved observations of the 
nearby PSQ J0330-0532 by Sanmartim et al.; ( 2014| >, although 
this object is not in our catalogue as it hosts a rather small pop¬ 
ulation of intermediate-age stars. 


5.2. Nuclear star formation 

We have seen that the intermediate-age stellar populations 
observed in K+A galaxies are not in place in PSQ, where 
intermediate-age stars are both less abundant and significantly 
metal poorer. We are forced to conclude, therefore, that either 
PSQ are not the immediate progenitors of K+A galaxies, or that 
the bulk of the young stars in PSQ are still being formed. This 
requires that vast amounts of metal-enriched molecular gas be 
funnelled onto the central regions during the merger in order to 
fuel a huge nuclear starburst, which is still not visible in PSQ. 
Unfortunately, the SDSS spectra of our objects lack the S/N and 
resolution required to infer the metallicities of the nebular gas 
in PSQ and discern nuclear stellar features, but we do know that 
in K+As the little nebular gas that remains has well over solar 
metallicities. 

The star formation rates of these putative nuclear starbursts 
would need to reach 10-100 solar masses per year (~ 10 10 M o 
in less than a few lOOMyr) and be extremely compact, perhaps 
extending over not much more than a few kpc (Sanmartim et 
|ak| |2014[ ). Under these conditions it would seem reasonable to 
expect that much of the luminosity that we observe as an ob¬ 
scured power-law component in PSQ would be produced not by 
an AGN, but by a very massive nuclear starburst. At an age of 
a few Myr the bolometric mass-to-light ratio of a simple star- 
burst is about M/LboI ~ 3 x 10 4 ( Leitherer 1999), so the mean 
bolometric luminosity of the power-law component in our PSQ, 
LboI = 10 44 9 erg s -1 , would imply a total mass of about 10 8 M a 
in young stars, about two orders of magnitude lower than the 
’’missing” intermediate-age stellar population in PSQ. 

An additional challenge for this scenario is that the young 
stars must form over time scales much shorter than lOMyr, 
which requires a very compact concentration of molecular gas. 
Compact nuclear starburst clusters may resemble type-I AGN as 
a consequence of supernovae evolving in a high density inter¬ 
stellar medium (Terlevi ch et al.| [l992), so a critical test of this 
Starburst-AGN connection is that at least some PSQ should have 
huge concentrations of cold, metal-rich molecular gas in their 
nuclear regions making them extremely luminous CO sources. 


5.3. Comparison with previous investigations 


Canalizo & Stockton (2001) have studied a sample of transi¬ 
tional AGN which seem to be objects in transition between the 
(U)LIRG and the PSQ phase. In Canalizo & Stockton ( {2013} ) 
they extended their study to more normal QSOs and consistently 
find that the intermediate-age populations of their objects - that 
they associate with remnants of major mergers - are offset from 
the nuclei by typically ~ 8 kpc. This is consistent with our con¬ 
jecture that the intermediate-age populations in PSQ must be the 
remnants of starburst clusters in the arms of the merging galax¬ 
ies. 


Our results are in broad agreement with the study of Cales 
|et al.| ( |20T3} who observed 38 PSQ selected from the SDSS, in 
particular with our conclusion that the putative AGN in PSQ co¬ 
exist with massive bursts of star formation. Unfortunately, there 
are only 12 objects in common between our two samples. Ten of 
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the objects in the Cales et al.]( |2013) sample are not in the DR7Q, 
while their selection criterion EW(Hh) >1A is substantially less 
stringent than our criterion of EW( Hdj >3 A, so it admits objects 
that we would not classify as PSQ. 

While as expected the majority of our objects fall in the 
strong AGN (sAGN) domain of the WHAN diagnostic diagram, 
most of the objects in the Cales et al. sample appear to be ei¬ 
ther LINERS or star-forming, with the caveat that they use the 
[OIII]/H/J vs. [Nil]/Ha BPT plane to classify their objects. In 
PSQ H/J is rather weak and extremely sensitive to corrections 
for underlying absorption, so using BPT diagrams involving H J3 


seems risky. However, the BH masses derived by Cales et al. 
( j2013| > using a completely different methodology are substan¬ 
tially similar to the BH masses of our PSQ, and this indicates 
that, even thought the stellar populations may differ, the AGNs 
themselves seem to be similar. 

In their study of J0330-0532, Sanmartim et al. (2014) found 
that the intermediate-age population avoids the nuclear regions 
and is distributed along a disk-like structure, which would be 
consistent with our results, although as mentioned above, this 
object does not meet our selection criterion for PSQ. 


6. Conclusions 

We have compiled a sample of 72 bona fide post-starburst QSOs 
(PSQ) selected from the SDSS DR7 quasar catalogue purely on 
the basis of the equivalent width of the Balmer H<) lines. We have 
used published photometry and SDSS spectroscopy to study the 
stellar populations, gas properties, and SEDs of the sample. Our 
main conclusion is that PSQ do not (yet) contain the large pop¬ 
ulation of metal-rich intermediate-age stars that we observe in 
post-starburst (K+A) galaxies. Thus, either PSQ are not the im¬ 
mediate progenitors of K+As, or they contain hugely massive 
nuclear starbursts where the bulk of their intermediate-age stars 
are still being formed. 

Our results, however, are based on rather low S/N SDSS 
spectra of these relatively faint and distant objects. Further 
progress in the understanding of PSQ, therefore, requires sub¬ 
stantially deeper observations in order to refine the population 
synthesis models and to investigate the kinematics and physical 
conditions of the nebular gas as well as to peek into the nuclear 
regions in search for the hidden AGN. 

The majority of our PSQ are not strong radio or X-ray 
sources, although the number of objects observed in X-rays re¬ 
mains rather small. Thus, deep X-ray and mm wave observa¬ 
tions are critically needed to establish the nature of the AGN in 
these objects and to check for the presence of nuclear starbursts. 
The present paper reports the initial steps of such research pro¬ 
gramme. 
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